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ABSTRACT 


Reservoir fluid properties are very important in reservoir engineering computations such as material balance calculation, well test 
analysis, reserves estimate and numerical reservoir simulations. Ideally, these properties should be obtained from laboratory pressure- 
volume-temperature (PVT) analysis. Quite often, however, these measurements are either not available, or very costly to obtain. In such 
cases, empirically derived correlations are used to estimate the needed properties, all computation therefore, will depend on the accuracy 
of the correlations used for estimating the fluid properties. Hence in this study, Standing's correlation for estimating the solution gas- 
oil ratio was optimized using a particle swarm optimization (PSO) algorithm to minimize the error associated in estimating solution 
gas-oil ratio from correlation at various depletion pressure. The optimized correlation was taken as a function of bubble point pressure, 
API gravity, gas gravity and reservoir temperature. PVT data from differential liberation test was used to validate this study's correlation 
and the result obtained shows that the optimized correlation for this study matches closely with the experimental values, also the newly 
optimized correlation was validated with other models and the results gave the least average relative error of 3.34 and a correlation 


coefficient of 0.998 after 216th successive iterations by the particle swarm optimization algorithm. 
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INTRODUCTION 

Pressure-volume-temperature (PVT) analysis is the study of the behavior of vapor and liquid in petroleum reservoirs 
in terms of phase behavior and composition. PVT properties are used to estimate reserves, evaluate the performance 
of oil and gas reservoirs, design production facilities and production operations. Also, to estimate stock tank gas 
production rate, cumulative gas-oil ratio, stock tank gas-oil ratio is a key input parameter (Hillary and Okotie, 2016). 

Determination of PVT properties among other data, for estimation of the stock oil initial in place and the 
evaluation of the future performance of hydrocarbon reservoir, is usually associated with uncertainties and as such, it 
is vital to properly characterize these properties which are key input data for most oil and gas calculations. According 
to Ikiensikimama and Egbe (2006), at the earlier stages of a well, it can be difficult or economically impractical to 
obtain reliable measurements of PVT data. In a scenario where the fluid samples are available, they can be subjected 
to PVT analysis to determine their properties, but samples are often suspected and PVT analysis usually applies only 
at reservoir temperature. In addition, Ikiensikimama (2008) has it that the reserves estimation and the design of the 
best depletion strategies are only feasible when a realistic and precise values of reservoir fluid properties are available. 

As stated by Okotie et al. (2017), to estimate the PVT properties of a reservoir fluid, the fluid is usually 
sampled and taken to the laboratory for experimental analysis such as saturation pressure (Dew point) at reservoir 
temperature, constant composition expansion test for black oil and compositional reservoir fluid to determine 
properties such as relative volume, vapor z-factor and liquid drop out. Differential Liberation/Vaporization test for 
black oil only to determine vapour Z factor, liquid density, gas-oil ratio, relative volume (formation volume factor), 
gas gravity, liquid viscosity, vapor viscosity. constant volume depletion test for compositional reservoir fluid only to 
determine the retrograde liquid drop out, cumulative fluid produced, vapor z factor, specific gravity of produced fluid 
plus, mole weight of produced fluid plus, final weight of produced fluid plus, produced vapor composition and finally, 
separator test for black oil and compositional reservoir fluid to estimate the gas-oil ratio and stock tank formation 
volume factor. 

Carrying out this study in the laboratory is usually expensive and time-consuming, thus, Engineers in the 
field resorted to using existing correlation to estimate these properties. The results from these existing correlations are 
approximations of the field data which yield a considerable amount of error, in order to minimize error associated in 
estimating the solution gas-oil ratio from experimental, an optimization algorithm is required. Thus, this study is aimed 
at minimizing the error between the experimented results and the result from existing correlations for accurate 


estimation of solution gas-oil ratio from correlations 


EXISTING PVT CORRELATIONS 
Standing (1947) correlation for gas-oil ratio (GOR) is given by equation (1): 
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D x 1.2048 
R, = y; | nat 1.4) 10 ] (1) 
X — 0.0125API — 0.00091(T — 460) (2) 
Oloruntoba and Onyekonwu GOR equation is given by equation (3): 
R, ane 1070:59255 Qo ae pert (3) 


Glaso (1980) proposed a correlation for estimating the gas solubility as a function of the API gravity, pressure, 


temperature, and gas specific gravity given by equations (4-5): 


AP]°-289 1.2255 
R, — yg (a = oad 10°] (4) 


X = 2.8869 — (14.1811 — 3.3053 log P)? (5) 





Al- Marhoun (1985) developed a correlation for estimating the solution gas oil ratio for Middle East crude 
oil. These correlations were based on a database of 69 bottom hole fluid samples and expressed as a function of 
reservoir temperature, gas gravity and oil gravity given by equation 6: 


1 
` 1.87784 
PYg 0.715082 


R: = | 9,00538088y,24437(T + 459,67)132657 (6) 


Petrosky and Farshad (1993) used a nonlinear multiple regression software to develop a gas solubility correlation. The 
authors constructed a PVT database from 81 laboratory analyses from the Gulf of Mexico crude oil system. Petrosky 


and Farshad proposed the following equations (7-8): 


1.73184 


fi [G5 + 12.340) y29*2910*] (7) 


X = 7.916(10-*) (API):5*1? — 4.561(10-5)(T — 460)13?11 (8) 


Kartoatmodjo and Schmidt (1994) developed a new set of empirical correlation based on a large collection of data 


developed from all over the world. Their equations (9-10) are given by 


R, c 0.0599 org o p aaa, iin rime 1) (API ae 30 (9) 
R, = 0.03150y 27997 pro LO eer are T4898) (API S. 30) (10) 


16 | Okotie and Fasanya, 2020 


0 0 
ij» i 


Initialize position x;;, C1,C2, velocity Vip 


0 
evaluate fj; 
D=max number of dimension, P=max number of particles, N=max number of iterations 


ti) 
Choose "jT; 





t+1 _ t t t t t t 
Vg HV + Gain [Piest i — x] + Cori; Goose = xi] 






xs = x + yt 


J«D 


NO 


Evaluate fj using xj; 


Yes 
t-— t ub 
fü — foesti Dbesti = Xij 
NO 


Yes 


t 
fij = fabest A fapest = fi Pbesti x 








NO 


un 


NO 


Figure 1: workflow algorithm for the PSO (gbest PSO) 
(Source: Satyobroto et al. (2001)) 


17 | Okotie and Fasanya, 2020 


THEORECTICAL CONCEPT OF THE PSO ALGORITHM 
Step 1: Choose the number of particles 

Step 2: Initialize the initial positions of the particles 

Step 3: Evaluate the objective function at the initial positions 
Step 4: Set the iteration number as t = i1 

Step 5: Find the personal best for each particle 

Step 6: Find the global best 

Step 7: find the velocities of the particles 

Step 8: Find the new values of the particles position 

Step 9: Find the objective function values of step 6 


Step 10: Stopping criterion: 


If the terminal rule is satisfied, go to step 4, otherwise stop the iteration and output the results. 











INPUT DATA 
Table 1: Differential liberation test data 
SAMPLE No 1- SEPARATOR TEST 
(T=129°F, P,=2405psig) 
Rysp (cuft @ 60°F and 14.65psia per STB @ 60°F) 701 
Bosp (barrels of saturated oil @ P, and T per STO @ 60°F) 1.356 
Stock Tank Gravity (API @ 60°F) 37 
Yg (from separator air = 1) 0.743 
Pressure(psig) Gas solubility (scf/STB) 
2405 737 
2200 684 
1950 620 
1700 555 
1450 492 
1200 429 
950 365 
700 301 
450 235 
200 155 
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RESULTS 
ESTIMATED RESULT FROM EXISTING CORRELATIONS 


Table 2: Result of gas solubility (scf/stb) from existing correlations 





Pressure(psig) Experimental Glaso Standing Petrosky and Fashad 





2405 737 665 703 685 
2200 684 500 735 508 
1950 620 425 658 467 
1700 555 388 575 398 
1450 492 319 493 342 
1200 429 285 309 297 
950 365 205 242 255 
700 301 165 185 200 
450 235 100 100 170 
200 155 55 45 125 





The results are plotted in Figure 2 to show the disparity in value of the correlations from the experimented value. 
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Figure 2: Result of gas solubility from correlations and experimental data 


ESTIMATED RESULT FROM PSO ALGORITHM 

The result of Standing's correlation gave the largest disparity than the other two correlations, hence, it was optimized 

with the particle swamp optimization algorithm to achieve the objective of this study. The particle swamp optimization 

algorithm is an iterative process that cannot be done manually, thus, it was programmed in Microsoft Excel to get the 

global best value for the constant A, B and C respectively as shown in Figures A1, A2 and A3 in the appendix section. 
Convergence criteria is a phenomenon of PSO in which all particles tend to converge to a single value as 


shown by Figures Al, A2 and A3 respectively where all three particles maintain a single value at the 216" iteration 


19 | Okotie and Fasanya, 2020 


number. These single values represent the optimum solution to the objective function as discussed above; the values 
of x; = (A; Bi Ci) at the 216" iterations are x216 = (17.53468968, 0.053651478, 0.745392522) respectively. 
Hence, the new empirical correlation used in estimating the gas solubility derived from particle swam optimization 


(PSO) algorithm, an optimization of standing correlation is given by equation | & 2. 


p 0.745392522 
+ 1.4) 10 a reponat ede] (14) 


Rsi optimized = Yg (3753468568 


Table 3 shows the estimated result obtained from the new correlation developed in this study. 


Table 3: Gas solubility (scf/stb) result for this study 











Pressure (psig) Experimental Glaso Standing Petrosky and Fashad (This study) 
2405 737 665 703 685 724 
2200 684 500 735 508 678 
1950 620 425 658 467 620 
1700 555 388 575 398 560 
1450 492 319 493 342 498 
1200 429 285 309 297 434 
950 365 205 242 255 365 
700 301 165 185 200 294 
450 235 100 100 170 214 
200 155 55 45 125 122 
RESULT OF CROSS PLOT 


To graphically illustrate the success of the existing correlations, a performance or cross plot was performed as shown 
in Figures 4 — 6, which is a plot of the estimated values versus the experimental values. A 45° straight line was drawn 
on the cross plot on which the estimated value is equal to the experimental value. Fitting the straight line at 45° clearly 
indicates that, the closer the plotted data points are to this line, the better the correlation. The result of Petrosky-Fashad, 


Glaso, and Standing correlations show a clear disparity from the experimented. 
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Figure 3: Cross plots for Gas Solubility using Standing Figure 4: Cross plots for Gas Solubility using Glaso 
Correlation Correlation 
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Figure 5: Cross plots for Gas Solubility using Petrosky Figure 6: Cross plots for gas solubility for this study 


and Farshad correlation 


STATISTICAL ERROR ANALYSIS 

According to the Tables 4, it is clarified that the errors calculated by this study’s correlation are lower than the other 
three correlations as shown in Figures 3-6 and newly developed correlation is most accurate for pressure below the 
bubble point pressure, yielding an average percent relative error of 3.34% and correlation coefficient of 0.998 which 
means that the sum of squares about regression is minimal i.e. the unexplained variation by regression and thus agree 


that the best correlation should have the least average percent relative error and highest correlation coefficient between 


zero and one. 
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Table 4: Statistical analysis of correlations 











Statistical parameter This study Standing Glaso Petrosky and Farshad 
correlation correlation correlation 

Average percent relative error 3.34 30.22 40.11 29.3 

(%) 

Sum squared residual 364985 566488 522387 373945 

Correlation coefficient 0.998 0.9167 0.8097 0.8089 

CONCLUSIONS 


Based on the results of this evaluations, the following remarks are worth mentioning: 


An empirical correlation for estimating the gas solubility at the bubble point pressure and below for black 
oils has been developed. 

The result obtained with the data from the differential liberation test, shows that the correlation developed 
in this study performs better than Standing’s, Glaso’s and Petrosky’s correlations. Hence, engineers can 
rely on the newly developed correlation to an extent after the have validated it with their field data 
because it is stated in literatures that correlation performs better in the region it was developed. 

The statistical result indicates a lower value of average relative error and a better coefficient of correlation 
for this study than for Standing’s, Glaso’s and Petrosky’s correlations at pressures below bubble point. 
The particle swarm optimization tool achieves a better accuracy by minimizing the objective function 
generated in this study. Above bubble point, solution gas oil ratio is the same with solution gas oil ratio 
at bubble point. 

Further work is recommended in developing a modification of Particle swarm optimization with faster 


convergence rate. 
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APPENDIX 
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Figure A1: Constant A values from PSO iteration 
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Figure A2: Constant B values from PSO iteration 
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Figure A3: Constant C values from PSO iteration 
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